Magnetorotational-type Instability in Couette— Taylor Flow 
of a Viscoelastic Polymer Liquid 
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We describe an instability of viscoelastic Couette-Taylor flow that is directly analogous to the mag- 
netorotational instability (MRI) in astrophysical magnetohydrodynamics, with polymer molecules 
playing the role of magnetic field lines. By determining the conditions required for the onset of 
instability and the properties of the preferred modes, we distinguish it from the centrifugal and elas- 
tic instabilities studied previously. Experimental demonstration and investigation should be much 
easier for the viscoelastic instability than for the MRI in a liquid metal. The analogy holds with 
the case of a predominantly toroidal magnetic field such as is expected in an accretion disk and it 
may be possible to access a turbulent regime in which many modes are unstable. 
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The magnetorotational instability (MRI) is widely re- 
garded as one of the most important instabilities in as- 
trophysics. It occurs in a differentially rotating fluid of 
sufficiently high electrical conductivity in the presence 
of a magnetic field, when the angular velocity decreases 
with distance from the axis of rotation. The MRI is espe- 
cially relevant to astrophysical accretion disks, where it 
is believed to give rise to turbulence and angular momen- 
tum transport, but it can also occur in stellar interiors 
[3, 0, and references therein] . 

The mechanism of instability can be visualized by con- 
sidering two fluid elements connected by magnetic field 
lines, which act as an elastic tether. The elements ini- 
tially orbit at the same radius and one is displaced inward 
so that its angular velocity increases. Being retarded by 
the tether, it loses angular momentum and moves to a 
still smaller orbit so that the process runs away. Both 
inertial forces and high conductivity are important to 
the mechanism. 

Experimental demonstration of the MRI is regarded as 
highly desirable but has proved difficult because of the 
need to achieve a sufficiently large magnetic Reynolds 
number (Rm = LU^lqct, where L and U are charac- 
teristic length and velocity scales of the flow and a is 
the electrical conductivity). Even for the most highly 
conducting liquid metals, the required Reynolds number 
(Re = LU/v, where v is the kinematic viscosity) is ex- 
tremely large, in excess of 10 6 . 

A natural choice of experiment is the Couette-Taylor 
flow of a liquid metal between differentially rotating con- 
centric cylinders. By setting the angular velocities of the 
cylinders in a "Keplerian" ratio, Q oc r~ 3 / 2 , and intro- 
ducing a uniform axial magnetic field, the MRI becomes 
possible while the centrifugal (inertial) instability is sup- 
pressed according to Rayleigh's criterion The water 
prototype 0] of a planned experiment using liquid gal- 
lium [f| has attempted to deal with the strong end-effects 
that arise at large Re when the cylinders are not much 



taller than the gap between them. 

Behavior related to the MRI may have already been 
observed in a spherical Couette flow of liquid sodium 
Q . The interpretation of this experiment is difficult be- 
cause of the presence of hydrodynamic turbulence and 
a complicated pattern of differential rotation in the un- 
magnetized state. A different instability, which relies on 
the presence of a helical magnetic field, has also been 
reported in a cylindrical experiment 0, S| ■ This result 
is controversial [9| and, since the instability occurs at 
Rm <C 1, it does not resemble the MRI as conventionally 
understood. 

In this Letter we describe an instability that is directly 
analogous to the MRI and occurs instead in a simple vis- 
coelastic liquid, with polymer molecules playing the role 
of magnetic field lines. We suggest that an experiment 
to study this MRI analog would be much easier and may 
in some respects represent the astrophysical systems of 
interest more closely than liquid metal experiments. 

The characteristic properties of magnetohydrodynam- 
ics (MHD) in the limit of large Rm are that the magnetic 
field lines are advected and distorted by the fluid flow and 
can be identified with material lines, and that the field 
exerts a tension through the Lorentz force. These two 
properties imply that the magnetic field imparts an elas- 
ticity to an electrically conducting fluid, which allows the 
propagation of Alfven waves, for example. 

Directly analogous behavior occurs in viscoelastic flu- 
ids such as dilute solutions of long-chain flexible poly- 
mer molecules in water or organic solvents. The poly- 
mer molecules also impart an elasticity to the fluid be- 
cause they are advected and distorted by the flow and 
exert a tension. In earlier work [l~oj | we drew attention 
to the physical and mathematical similarities between 
these situations. The mathematical viewpoint is that the 
force per unit volume is the divergence of a symmetric 
stress tensor which is Lie-transported by the fluid flow in 
the sense of a second-rank contravariant tensor field (its 
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upper-convected derivative is zero). This result applies 
asymptotically both to the (Maxwell) magnetic stress in 
MHD in the limit of large Rm and to the polymeric stress 
in the limit of large Deborah number (De = St, where 
S is the shear rate and r is the relaxation time of the 
fluid). 

We also showed [l(| that there is a direct analog of the 
MM in rotating shear flows of a viscoelastic fluid, using 
a Cartesian model that can be considered as a local ap- 
proximation of a Couette-Taylor flow when the rotation 
and shear rates are comparable. For De 3> 1 the poly- 
meric tension is parallel to the flow and corresponds to an 
equivalent magnetic field also in that (azimuthal) direc- 
tion. Growth rates and eigenfunctions of the viscoelastic 
and MHD problems were found to be very similar. 

Here we present a new analysis of the linear stabil- 
ity of viscoelastic Couette-Taylor flow, applicable to an 
experiment with cylindrical geometry. By optimizing 
the growth rate over all wavenumbers, we determine the 
likely initial outcome of such an experiment. We also 
show how the MRI analog can be distinguished from 
other types of instability in the same system. 

A viscoelastic fluid is governed by the equation of mo- 
tion, 



' du _ 



= -Wp + V - T + p s V 2 u, (1) 



the incompressibility condition, V • u — 0, and, in the 
case of an Oldroyd-B fluid, the constitutive equation 



u • VT — ( Vu) T • T — T • Vu 



"<9T 

= n P [Vu + ( Vu) T ] 



(2) 



Here p, p, and u are the density, pressure, and velocity of 
the fluid, fi s and p, p are the solvent and polymer viscosi- 
ties, and r is the relaxation time of the polymeric stress 
T. The widely adopted Oldroyd-B model is the simplest 
nonlinear description of a dilute polymer solution and 
can be derived directly from a kinetic model [ll[ . 

In the Couette-Taylor experiment the basic flow is u = 
rf2(r) e<jj, where (r, cj>, z) are cylindrical polar coordinates, 
fi(r) = A + Br~ 2 , and the constants A and B are such as 
to match the angular velocities Q\,2 of the inner and outer 
cylinders of radii n,2- The polymeric stress components 
are T r( f, — ~2p p Br~ 2 and = 8p p TB 2 r~ 4 , the latter 
being analogous to the "hoop stress" of a magnetic field 
B^> — (poT^) 1 / 2 . End-effects should be unimportant 
provided that the height of the cylinders is much larger 
than the gap width d — — r± . 

Two dimensionless numbers associated with the exper- 
iment are the Reynolds number Re = r\d\U,2 — Qx\/v, 
which measures the relative importance of inertial and 
viscous effects (here v — (/z s +/x p ) /pis the total kinematic 
viscosity), and the Deborah number De = T"i|f22 — ^i\r/d, 
which quantifies the importance of elasticity. 



We linearize equations ([I]) and J2]) about the basic flow 
and seek solutions with the dependence exp(si + im<j) + 
ikz), obtaining a sixth-order system of ordinary differen- 
tial equations for the radial structure of an eigenmode 
with complex growth rate s, integer azimuthal mode 
number m and real axial wavenumber k. No-slip bound- 
ary conditions apply at r — r\^- 

We solve the problem by two independent methods. 
First, we apply Chebyshev collocation on a Gauss- 
Lobatto grid to convert the system into an algebraic gen- 
eralized eigenvalue problem, which we solve directly using 
a LAPACK routine. Second, we use a shooting method 
based on a 4th/5th order Runge-Kutta integrator with 
adaptive stepsize. The two routines agree to very high 
precision. The Chebyshev method is useful for examining 
all the eigenmodes with a sufficiently simple radial struc- 
ture to be adequately resolved (spurious modes are also 
generated which must be rejected). The shooting method 
is good for following a single mode as the parameters are 
varied. 

At each point in the (De, Re) plane we seek to maxi- 
mize the growth rate Re(s) with respect to both m and 
k. Once the preferred wavenumbers are obtained at one 
point, the optimal solution is tracked through the plane, 
making the necessary adjustments to m and k. 

Two different instabilities have previously been re- 
ported in viscoelastic Couette-Taylor flow 12fl . One is 
the centrifugal (inertial) instability (CI) of Rayleigh and 
Taylor which occurs at sufficiently large Re even in 
the Newtonian limit (De = 0). The other is the purely 
elastic instability (EI) of Muller at al. [l3} and Larson 
et al. 14|, which occurs at sufficiently large De even in 
the inertialess limit (Re = 0). The CI is driven by shear 
energy and requires the specific angular momentum \r 2 fl\ 
to decrease with r, while the EI is driven by elastic energy 
and simply requires shear and curved streamlines. 

The MRI analog, however, depends on both elasticity 
and inertial forces, and therefore requires both Re and 
De to be sufficiently large. It is driven by shear energy 
and requires the angular velocity 101 to decrease with r 
1- 

We tested our methods by reproducing standard re- 
sults on the CI of a Newtonian fluid and also the results 
of Ref. [l(| in the limit of a narrow gap. 

We made strenuous efforts to reproduce the results of 
Ref. (HI, where a problem identical to ours was consid- 
ered in a different parameter regime, but we obtained 
only qualitative agreement with many of their findings. 
To verify our own results we derived the linearized equa- 
tions by two independent methods and also obtained 
matching numerical results with completely independent 
codes. 

For the purposes of illustration we adopt the param- 
eters T\/r% = 0.95 and p s /Pp = 1- The latter ratio 
could be either small or large in practice, depending on 
the choice of solvent. Boger fluids, which avoid shear 
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FIG. 1: Contours of the growth rate Re(s) /fii , optimized with 
respect to m and k, for case 1 (p = —3/2). Contour values 
increase from to 0.35 in steps of 0.025 starting from the 
thick line (negative contours are not plotted). The preferred 
azimuthal mode number m is indicated by numerals. 

thinning and conform more accurately to the Oldroyd-B 
model, have more viscous solvents. 

We set the angular velocities of the cylinders in the 
ratio ^2/^1 = ( r 2Ai) p , with p = —3/2 (case 1), p = 
+3/2 (case 2) or p = —3 (case 3). Case 1 (a Keplerian 
ratio) should exhibit both EI and MRI, case 2 only EI, 
and case 3 all three instabilities. 

In Fig. 1 we show the optimized growth rate in the 
(De,Re) plane for case 1. The EI is seen at Re = for 
De > 19; it prefers a small value of m and its growth rate 
is limited here to < 0.1 Oi. As Re is increased, instabil- 
ity occurs for a much wider range of De and with larger 
growth rates. This is the MRI analog. The preferred 
value of m increases with Re/De as the equivalent mag- 
netic field becomes weaker [lCj. The EI prefers a larger 
value of k. As expected, the CI is absent and there is no 
instability in the Newtonian limit De = 0. 

A typical eigenfunction of the MRI analog at the on- 
set of instability is shown in Fig. 2. This solution is 
qualitatively similar to Taylor vortices and has a vertical 
wavelength comparable to the gap width but is nonax- 
isymmetric and would appear as a rotating spiral pattern. 
However, nonlinear pattern selection may favor a linear 
combination of upward and downward spirals, resulting 
in a "ribbon" structure. 

In case 2 (Fig. 3) the ratio of shear to rotation is similar 
to case 1 but of opposite sign, so that the MRI and CI 
are absent. We find that the EI operates in a very similar 
fashion at small Re but is suppressed by increasing Re 
at fixed De. A comparison of Figs 1 and 3 makes it 
clear that only the very low-Re behavior in Fig. 1 can be 
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FIG. 2: Eigenfunction of the preferred mode in case 1 at the 
onset of instability for (De, Re) = (3.3,40), where (m,k) = 
(11,95). The azimuthal velocity perturbation is indicated by 
color/grayscale. One vertical wavelength is shown. 




FIG. 3: Growth rate for case 2 (p = +3/2). Contour values 
increase from to 0.07 in steps of 0.005. The preferred value 
of m is 2 throughout. Note the different horizontal scale. 

attributed to the EI. 

In case 3 (Fig. 4) we observe similar behavior to case 1 
in much of the (De, Re) plane. However, the CI is seen 
operating at small De for Re > 85, where it prefers m = 
(the classic axisymmetric Taylor vortices). Fig. 4 there- 
fore illustrates EI, MRI, and CI, but most of the un- 
stable region can be identified as MRI by comparison 
with Figs 1 and 3. The growth rates and preferred mode 
numbers of the MRI are somewhat different from case 1 
because of the different ratio of shear to rotation. Wc 
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FIG. 4: Growth rate for case 3 (p = —3). Contour values 
increase from to 0.7 in steps of 0.05. Note the different 
vertical scale. 



manifestation of the MRI analog. 

The asymptotic nature of the analogy between vis- 
coelastic and MHD flows [l(| means that it is more accu- 
rate for larger values of Re and De. At the onset of insta- 
bility, as in Fig. 2, the analogy is good but not perfect, 
ft is better in the upper right part of Fig. 1, where many 
modes are unstable and the flow may become turbulent. 
This situation would be analogous to a turbulent accre- 
tion disk containing a predominantly toroidal magnetic 
field. In contrast, liquid metal MRI experiments employ 
an axial magnetic field and are not expected to be able 
to access a regime of MHD turbulence. In these respects 
the polymer experiment may offer a more faithful repre- 
sentation of the astrophysical systems of interest. Other 
major advantages are that the materials are cheap and 
safe, and that optical visualization or velocimetry can be 
carried out. 

This work was supported in part by STFC. 



find that the three instabilities are described by a mode 
that changes its character in a continuous way across the 
parameter space. 

Experimental and theoretical studies of instabilities in 
viscoelastic Couette-Taylor flow have been carried out 
since the 1960s but have not revealed the three types of 
characteristic behavior seen here. The MRI analog is ex- 
hibited most clearly, as in our case 1, when the angular 
velocity decreases outward while the angular momentum 
increases. The few experiments in which the inner and 
outer cylinders can rotate independently (e.g., [Hj]) ap- 
pear not to have examined this regime. 

Previous studies have described either a reduction of 



the critical Re for the CI as De is increased [12L and ref- 
erences therein] or a reduction of the critical De for the 
EI as Re is increased 17[ . The physical nature of this in- 
termediate "inertioelastic mode" has not been discussed. 
A similar behavior can be seen in case 3 (Fig. 4) but 
we interpret it as the MRI analog. The MRI tends to 
dominate the CI when both the angular velocity and the 
angular momentum decrease outward, unless De is small. 
In case 2 (Fig. 3), where the MRI is absent, increasing 
Re tends instead to suppress the EI. 

An experimental study of the MRI analog would re- 
quire a standard Couette-Taylor apparatus and a suit- 
able polymer solution. The polymer viscosity and relax- 
ation time can be adjusted by varying the concentration 
of the solution. Care must be taken to maintain the tem- 
perature of the fluid in the presence of viscous dissipation. 
Some effects such as shear-thinning in lower- viscosity sol- 
vents may require the use of a more sophisticated consti- 
tutive model than Oldroyd-B. Subcritical, hysteretic be- 
havior has been reported in viscoelastic Couette-Taylor 
flow [It], EH and this may also affect the experimental 
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